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Abstract 

We study a few problems related to Markov processes of flipping triangulations of the 
sphere. We show that these processes are ergodic and mixing, but find a natural example 
which does not satisfy detailed balance. In this example, the expected distribution of 
the degrees of the nodes seems to follow the power law 

1 Introduction 

We consider a Markov chain on triangulations of the sphere (or other surface). Let T 
denote the set of triangulations, by this we mean the set of all combinatorially distinct 
rooted simplicial 3-polytopes. 

Tutte [7J showed that their number is asymptotically 



as the number n of vertices goes to oo. Of course, Euler's theorem holds for such 
triangulations, and this means that when there are n nodes, there are also Srt — 6 links 
and 2n — 4 triangles. 

For an element T E T,we denote by M{T) the set of nodes and by C{T) the set of 
links. 

For any link I (connecting the nodes A and B), we consider the "complementary" 
link £' , which is defined as follows: if (A, B, C) and (A, B, D) are the two triangles 
sharing the link I, then I' is the link connecting C and D. 

We assume that for any r G a probability Pt^ is given on ZI(r), /.f?., P7^(£) — 
1. We define a Markov chain on T as follows. We first choose a link I G C(T) at 
random (with probability Vxit))- 

• If the link belongs to £(T), we do not change T and proceed with the next 
independent choice of a link. 

• If I' does not belong to C{T), we erase £ and replace it by We obtain in this 
way a new triangulation T' and we proceed with the next independent choice of a link. 
This replacement of I by is commonly called a. flip see [6], or a Gross-Varsted move 
0. See Fig. [1] 

We will denote by P(r'|T) the transition probability of this Markov chain. 




(1.1) 
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Figure 1: A flip: the link (A-B) is exchanged with (C-D). 

2 Properties of the Markov chain 

We now fix n and let Tn denote those triangulations with n nodes. 

Proposition 2.1 Assume that iidTeTn ^T^^eeccT) Pt(^) > 0. Then the Markov chain 
defined in Sec.^is irreducible and aperiodic. 

Proof. It is well known (see |6|) that by flipping links as described above one can 
connect any two triangulations of T„ (one shows that any T can be flipped a finite 
number of times to reach a "Christmas tree" configuration). 



n 




Figure 2: The "christmas tree" with n nodes, the "branches" between 5 and n not being 
shown. Any triangulation can be brought to this form by a sequence of flips. 

Since by our hypothesis any such (finite) succession of moves has a non-zero prob- 
ability this shows the irreducibility of the chain. To prove aperiodicity, we have to 
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prove that for any high enough iterate of the transition matrix, all the diagonal entries 
are positive. By the previously mentioned result, it is enough to show that we can con- 
struct cycles of length two and three for the "Christmas tree". Cycles of length two are 
easily obtained by flipping a link back and forth. For cycles of length three, we con- 
sider the sub "Christmas tree" of size six at the base of the complete "Christmas tree", 
see Fig.|2] We enumerate the nodes as in the figure, assuming n > 7. In particular node 
3 has degree 3, nodes 4 to n have degree 4 and nodes 1 and 2 degree n — 2. 




Figure 3: The three stages in the cycle of 3 flips regenerating the Christmas tree with n 
nodes. 

The cycle of length 3 is obtained by performing the following flips 

(1-4) ^ (3-5) 
(2-3) ^ (1-4) 
(4-5) ^ (2-3) 

after which we get again the "Christmas tree" with nodes 3 and 4 exchanged. □ 

Remark 2.2 Of course, the condition of Prop. 12.11 is not necessary, but we do not 
know any simple other criterion in terms of the Pt, but one can think for example of 
conditions involving two successive flips. 

From this result we conclude that there is only one invariant probability measure, 
and with this measure the chain is ergodic and mixing. 



3 Two Examples 

The easiest example is that where one chooses a link uniformly at random. Then one 
gets the uniform distribution on T, and, using this simple fact, many properties of this 
process can be deduced, see, e.g., f^. 

Here, we consider another example, which was suggested to us by Magnasco [3^ 4iJ. 
This process consists in first choosing a node uniformly and then to choose uniformly 
a link from this node. Let n be the number of nodes. An easy computation, shown 
below, leads to 

where di{£\T) and d2{£\T) are the degrees of the nodes at the ends of link £ in the 
triangulation T. 

Proof. If ^ is a link, we denote by di the two nodes it connects. If £ is a link and i is a 
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node, we say that £ ^ i if i ^ d£. If i is a node, we denote by di(T) its degree in the 
triangulation T. We have from Bayes' formula 

PtW = I «)P(«) • 

Moreover, P(i) = 1/n for any i, Pxii | i) = if £ 7^ i and otherwise 
Therefore 

which is formula i3.1\ . □ 
It also follows directly from this expression that for any T G T, 



In this computation we have not used the fact that T is a triangulation. Therefore this 
relation holds for any graph. 

For the second model, we have 

Theorem 3.1 The Markov chain P( • | • ) w not reversible {when n > 7). 

Remark. We have not checked what happens for smaller n. 

In other words, one cannot easily guess the invariant measure from the transition 
probabilities. 

Proof. Assume the chain is reversible with respect to some probability P on T, namely 
for any T and T' in T we have 

P(T' I T)P(T) = P(T I T')P(r') . (3.2) 

If Ti , . . . , Tfc , Tk+i — Ti is any cycle of admissible flips, we must have 

A p(r,|r,+i) ^ 

lip(T,+i|T,) ■ 

We are going to show that there is a cycle of length 4 for the Christmas graph for which 
this is not true, see Fig.|4] 

Consider the following cycle for the "Christmas tree" with the same notations as 
before 

(1-4) ^ (3-5) 
(2-5) ^ (4-6) 
(3-4) ^ (2-5) 
(5-6) -> (1-4) 

An easy computation leads to 

Yt P(Tj\T,+i) ^ 10 

lip(T,+i|r,) 9 ' 

if the number of nodes is larger than 6. □ 
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4 Numerical simulation 

We have performed extensive simulations on the model described above. In this section 
we summarize the numerical findings, but the reader should note that we have no theo- 
retical explanation for the results. The main insight is that the model with the uniform 
measure |2| leads to an exponential degree distribution, while the model of |4| leads to 
a power law distribution in a sense which we make clear now, see Fig|5l 
We formulate the results as 

Conjecture 4.1 There is a probability measure p on the integers larger than 2 such 
that the average number of nodes of degree d divided by n converges when n tends to 
infinity to p(d). Moreover p has polynomial decay in the sense that d^'^p(d) converges 
to a nonzero finite limit when d tends to infinity. 

Remark 4.2 It should be noted that several deviations from a pure power law are 
present in these experiments and will not go away with large n. First of aU, nodes 
of degree 3 are less frequent than would be suggested by a power law. We attribute 
this to the impossibility of doing a flip if a node of degree 3 is chosen: All its edges 
are unflippable. Second, if there are n nodes, assuming an approximate power law of 
N{d) = c • we find c « 50n (from c^^J^g — n). Thus there should be no nodes 
for which N{d) < 1, that is SOnd^"* < 1 or d > (50n)^/^. However, the experiments 
clearly show the presence of "outliers" of much larger degree. Closer analysis (with, 
e.g., logarithmic binning) reveals that these outliers are spaced at equilibrium in a way 
to continue the measured power law. 
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Figure 5: A log-log plot for triangulations of size n = 8194, 32770, and 524290, after 
about 10^'^ flips. The data are the cumulated sum d(i) of number of nodes with degree 
> i. The straight part is well fitted with a law of d(i) ~ cd^^, so that the degree 
distribution seems to be ^ d^*. The outliers are produced by the lacunarity of the data 
when the expected number of nodes of a given degree starts to be less than 1 . 
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We have done extensive checks for the correlations of degrees of neighboring nodes 
in our simulations. Such correlations have been theoretically explained in [2 1 for the 
case of the uniform choosing rule. These correlations are difficult to measure, but no 
decisive deviation from independence was found, except for some obvious topological 
rules. 

There is a feeling in the community of specialists, be they interested in random 
triangulations, or in 2-d gravity (the dynamic dual of our problem) that the "typical" 
triangulation should be "flat" (which means that each node should (wants to?) have 
6 links). To measure the effect of the tails of distribution of degrees, we use combi- 
natorial differential geometry, as advocated by Robin Forman (Q, who introduces a 
notion of "combinatorial Ricci curvature" which, in our case of triangulations reduces 
to — 5di. Extensive simulations show that this quantity seems to grow more or 

less monotonically as the process reaches the equilibrium state. Note that since 
does not depend on the triangulations, we are just measuring the sum of the squares of 
the degrees. 

Another observation, which holds with very high accuracy is that once a node has 
been chosen, at equilibrium, exactly 50% of all attempted flips are not possible, because 
the "other" link is already present. This means that a tetrahedron is placed on top of 
a triangle. Note that the study of such "vertex-insertions" is already present in Tutte's 
work 1 7 1 . 
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